Hydrogen atom in strong magnetic field revisited 
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We derive in a straightforward way the spectrum of a hydrogen atom in a strong 
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' The spectrum of hydrogen atom in a strong magnetic field was found long ago pQ, 

and is presented now in textbooks (see, e.g., [2] (§112, Problem 3)). The approach used 
therein is as follows. At first, using the wave functions of electron in a magnetic field, one 
constructs an effective potential for the motion along the field, and then the spectrum in 
this potential is found. We present here a somewhat different solution of the problem, 
with the results coinciding in fact with those of Q3I2]- We hope however that our approach, 
physically straightforward and transparent, is of some interest. At least, the present note 
may be considered as a sort of mini-review on the subject. 

Our starting point is the obvious observation that in a sufficiently strong magnetic 



field H (the exact criteria are discussed below) the motion of an atomic electron becomes 
almost one- dimensional, along the magnetic field, in the Coulomb potential —e 2 /z. 
The corresponding one-dimensional wave equation is 



4 z. 

We have introduced in it the usual dimensionless variable: 

2z 



av ' 



here a = h 2 /m e e 2 is the Bohr radius, m e is the electron mass, v is the effective quantum 
number, related to the electron energy as 



^ m e e 4 . , 

K = ^7 • 2 

Equation J1J coincides exactly with the radial equation for the s-wave in the three- 
dimensional Coulomb potential — e 2 /r, and has therefore the common hydrogen spectrum 

E '<-W^ "= 1 ' 2 < 3 "" (3 > 

and the set of solutions 

u~(z) = exp(-z/2) zF(-n, 2; z) , v = n, (4) 
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where F is the confluent hypergeometric function (here and below we are not interested in 
the normalization factors). These solutions vanish at the origin and are trivially continued 
to z < 0. Thus obtained solutions on the whole z axis are odd under z — > — z (as reflected 
by the superscripts "minus" in (jSJ), (HJ)). 

There is however an essential difference between the present problem and the s-wave 
Coulomb one. In the last case (jlj) is the only solution. The reason is well-known. Naively 
the radial wave equation for R(r)(= u(r)/r) has two independent solutions, which behave 
for r — > as R ~ const (u ~ r) and i? ~ 1/r (u ~ const), respectively. However, in 
fact i? ~ 1/r is no solution at all for the homogeneous wave equation if the point r = 
is included, since A (1/r) = — 47ro~(r). As to our problem, equation ((TJ does not describe 
really the vicinity of z = since therein we have to consider seriously the magnetic field 
itself. Therefore, there are no reasons to discard those solutions of (0) which tend to a 
constant for small z (and of course decrease exponentially for z — > oo). 

Such solutions are presented in a convenient form in [2 a (Mathematical Appendices, 
§ d, (d.17)). To our purpose they can be written for z > as 

u+(z) = exp (-z/2) 1 1 - vz (in z F(l - u, 2; z) 

« r(l- y + fc)Ml- y + fc)-^(fc + 2)-V>(A; + l)] fc \) 
+ t r(l- 1/) *!(*+!)! "J/' (5) 

here ^(a) denotes the logarithmic derivative of the gamma function: ip(a) = r'(a)/r(a). 
Being trivially continued to z < 0, thus obtained solutions on the whole z axis are even 
under z — > — z (as reflected by the superscript "plus" in (jSJ) and in the corresponding 
eigenvalues below). 

Under any reasonable regularization of the logarithmic singularity at z — > 0, the even 
solutions should have vanishing first derivative at the origin. In this way we obtain the 
following equation for the eigenvalues of v: 

In— = ^ +if>{l-v); (6) 



here a# = yhc/eH is the typical scale for the radius of electron orbits in the magnetic 
field H. We are working in the logarithmic approximation, i.e. assume that 

A = In a/a H >1- (7) 

This allows us to use a crude cut-off at an for the formal logarithmic divergence at z — » 0, 
as well as to simplify somewhat this equation. 
The smallest root of equation (0) is 



which gives the ground state energy 
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Other roots of equation (0) are 



n+\ n = 1,2,3..., (10) 
A 



with the corresponding energies 



2h 2 n 2 




Let us mention that the one- dimensional Coulomb problem was considered in [3 - 5] 
with various regularizations of the singularity at z — > 0, but without any relation to the 
problem of the hydrogen atom in a strong magnetic field. 

But let us come back to our problem. The resulting spectrum of the hydrogen atom 
in a strong magnetic field looks as follows. Each Landau level in this field serves as an 
upper limit to the sequence of discrete levels of the Coulomb problem in the z direction. 
This discrete spectrum consists of a singlet ground state with the energy given by formula 
(JHJ), and close doublets of odd and even states of energies given by formulae (j3J) and (fTTj) . 
There is also a continuous spectrum of the motion along z above each Landau level. 

This picture is valid for sufficiently low Landau levels, as long as the radius of a 
magnetic orbit is much less than the Bohr radius. Obviously, in a strong magnetic field 
this description fails for large magnetic quantum numbers, i.e., in the semiclassical region. 
Here we can estimate the orbit radius directly from the well-known spectrum of electron 
in a magnetic field (see, e.g., [2] (§112, Problem 1)): 

t eH , . . m + \m\ , , 

E = h iV + 1/2, N = n p + J— L; 12 

m e c 2 

here n p is the radial quantum number in the xy plane, and m is the angular momentum 
projection onto the z axis. Now the semiclassical estimate for the magnetic radius is 



a H {N) w ^—N = a H VN. 

Thus the present picture of levels holds as long as 

A = In— > In N. (13) 
a H 

At last, let us consider the correspondence between the obtained system of levels in 
a strong magnetic field and the hydrogen spectrum in a vanishing field. A beautiful 
solution of this problem was given in (and is quoted in [T]). We would like to present 
the solution here as well (having in mind in particular that our note can be considered as 
a mini-review). 

The crucial observation made in [S| is as follows. While changing the magnetic field 
from vanishingly small to a very strong one, the number of nodal surfaces of a given 
wave function remains the same. A hydrogen wave function (in zero magnetic field) with 
quantum numbers n, I, m has n r = n — / — 1 nodal spheres and I — \m\ nodal cones with 
a z axis. With the increase of the magnetic field, the nodal spheres become ellipsoids of 
rotation, more and more prolate, tending to cylinders in the limit of infinite field. The 
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correspondence between n r and n p (n p being the radial quantum number in the xy plane 
in a strong magnetic field) gets obvious from this picture: 

n p = n r . (14) 

The evolution of the hydrogen nodal cones is less obvious. However, due both to equation 
(jl4j) and to the conservation of the total number of nodal surfaces, the number n z of 
the nodes of a solution of equation ((*[)) should coincide with the number of nodes of the 
corresponding spherical function, 

n z = I — \m\. (15) 

In other words, I — \m\ nodal cones of a hydrogen wave function evolve into n z planes of 
constant z corresponding to the nodes of an eigenfunction of equation (|T4^) . 

And at last, let us note that during the whole evolution of the magnetic field, m 
remains constant. 

Let us consider now, for instance, the ground state in the magnetic field, with N = 
(see (H2J)). Obviously, it is degenerate, and its corresponding magnetic wave functions have 
n p = and m — 0, —1, —2 ... . Let us confine further to its lowest sublevel corresponding 
to the ground state solution of equation (JTJ), with n z = 0. According to the above 
arguments, the hydrogen ancestors of those wave functions should have n r = and 
I — \m\. In other words, these ancestors are: 

Is; 2p, m = —1; 3d, m = —2; and so on. 
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